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The current status of axion physics is presented. There still exists the axion 
window 10 9 GcV < F a < 10 12 GcV. The recent CAST solar axion search 
experiment on the axion-photon-photon coupling strength has to be improved 
by a factor of 100 to reach down to the region of superstring axions. The 
calculable B and m u = cases for strong CP solutions, and axino cosmology 
in SUSY extension of axion are also commented. 
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1. Introduction 

Modern cosmology needs dark matter and dark energy in the universe: 
f^CDM — 0.23, Ha — 0.73. There are several particle physics candidates 
for CDM: LSP, axion, axino, gravitino, LKP and other hypothetical heavy 
particles with some kind of Z2 symmetries. 

The old electroweak scale axion is the pseudo-Goldstone boson 1 arising 
from breaking the global Peccei-Quinn (PQ) symmetry. 2 The very light ax- 
ion is the invention from the need to solve the strong CP problem through 
PQ symmetry with electroweak singlet field(s). 3 ' 4 Superstring axions 13 ' 14 
may be in this very light axion category. The existence of instanton solu- 
tions in nonabclian gauge theories needs 9 vacuum, introducing a CP vio- 
lating interaction. 5 In the 9 vacuum, the physically meaningful interaction 
is parametrized by 9 

——F^F^ = {FF}; 9 = 9 qC D + 9 wcak 



where #qcd is the value determined from high energy scale and 



'weak 
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Arg. Det. M q is the one contributed when the electroweak CP violation is 
introduced. Here 9 is the final value taking into account the electroweak 
CP violation. For QCD to become a correct theory, this CP violation by 
9 must be sufficiently suppressed. A nonvanishing value 9 contributes to 
the neutron electric dipole moment d n . From the experimental limit, 6 we 
obtain the bound 

\d n \ < 0.63 x 10 25 ecm -> \9\ < KT 9 . 

Why is this so small? It is the strong CP problem. There are three types 
for the solution: (1) Calculable 6, (2) Massless up quark, and (3) Axion. 

One may argue that there were no strong CP problem in the beginning. 
In particular in 5D extension, since the instanton solution is the one in 
4D. I think this does not work or at best belongs to the calculable 9 type, 
because in the 4D effective theory one can always consider a 4D theory after 
integrating out the 5th coordinate. Let us briefly comment on two solutions 
first. 

• The Nelson-Barr type: 7 CP violation is introduced spontaneously. 
So, original Yukawa couplings are real. Spontaneous CP violation 
is introduced at high energy by introducing vectorlike heavy quarks 
so that they mix with light quarks. If the heavy vectorlike quarks 
are not introduced, the CP violation of light quarks originated by 
the high energy scale CP violation will be tiny due to the de- 
coupling theorem. Not to be affected by the decoupling theorem 
and to guarantee a tree level Arg. Det. M q = 0, specific forms for 
Yukawa couplings are assumed: SU(2)xU(l) breaking real VEVs 
appear only between F — F Yukawas, and CP violating phases 
in the VEVs appear only in F — R Yukawas, where F are the SM 
fermions and R are the heavy fermions. If heavy vectorlike fermions 
are integrated out, the effective Yukawa coupling structure of the 
low energy sector is of the Kobayashi-Maskawa form. 

• Massless up quark: Suppose that we chiral-transform a quark, 

q -> e 475 V 

It is equivalent to changing 9^9 — 2a. Thus, if it is allowed 
to have such a symmetry then strong CP problem is not present. 
The massless quark case belongs here. This solution was known 
from the very beginning of the strong CP problem but was not 
taken seriously because the up quark seemed to be massive. 8 The 
problem is, "Is m u — allowed phenomenologically?" The famous 
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up/down quark mass ratio from chiral perturbation theory (cPT) 
calculation is — = I. But physics below 100 GeV is more involved. 



There is the determinental interaction of 't Hooft, pictorially shown 
as Below the electroweak scale quarks obtain mass. Suppose, the 



up quark is massless. Then, there is no strong CP problem. But 
chiral perturbation theory can be done with instanton generated 
up quark mass from the above 't Hooft interaction, m u = mdm s /A, 
where A is at the QCD scale. So it is the problem whether the 
instanton calculus really gives the desired magnitude, in which case 
9 is still unphysical. In the community, still there is a disagreement 
on this issue: Kaplan and Manohar (KM), and Choi belongs to 
the positive group, 9,10 and Leutwyler (L) belongs to the negative 
group. 11 CP even observables do not see m u . From the figure, for 
example, we have m u>e ff — mdm s / A, m^e// — + m„m s /A ~ 
m u m s /A, m s ^ e ff ~ m s . But CP odd observables see m u . Is Z = 
m u /md small? KM shows from the 2nd order cPT Z ~ 0.2, and 
they could not rule out the m u = case. Explicitly, cPT has the 
Lj parameter in the term Li(M^U — MU^) , where M = 3 x 3 mass 
matrix and C/ = 3x3 matrix for meson fields. KM shows 



where Lg, L§ are another parameters in the cPT. On the other hand 
L attempted to compute L7, using the QCD sum rule for the SU(3) 
singlet pseudoscalar rf dominance (similarly to the vector meson 
dominance) , 

L7 ~ £7,77' — (~2—4) x 10 _4 (Gasser — Leutwyler coefficients) (2) 

with a notable sign difference from (1). If (2) were true, the case 
m u = is ruled out. But Choi 10 argues that if rj gets mass from in- 
stanton calculus, which is the modern wisdom on the U(l) problem 



9 





(1) 
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resolution, he can change the sign of (2) to 

L 7 ~ (3 - 8) x 10~ 4 . (3) 
So, we can have the possibility of m u = 0. 

In recent years, lattice calculation has been performed toward this 

issue. 12 In 16 3 x 32 lattice calculation, they obtain: 2L S — L 5 ~ 10~ 4 

and m u /md = 0.484 ± 0.027. If true, m u = is ruled out. 

I consider that the problem on m u = is not completely settled 

yet, even though m u ^ seems to be the majority opinion of the 

community. 

• These show that the axion solution is the most compelling solution 
which is discussed in the subsequent section. 

2. Axion 

The axion potential is of the form where the vacuum is shown as a bullet. 



The vacuum stays there for a long time, and oscillates when the Hubble time 
(1/i?) is larger than the oscillation period (l/m a ): H < m a . This occurs 
when the temperature is about 1 GeV. Axion is directly related to 8. Its 
birth was from the PQ symmetry whose spontaneous breaking introduced 
a dynamical degree, a pseudo-Goldstone boson called axion. But "pseudo- 
Goldstone" nature is specific in axion in that axion is a pseudoscalar a 
without any potential except that arising from, 

— Ur— FF = —{FF}. (4) 
32^ 2 F a F a X 1 V 1 

This kind of nonrenormalizable term can arise in several ways. The first 
important scale is F a , defining the strength of nonrenormalizable interac- 
tion. It can arise from higher dimensional fundamental interactions with 
the Planck scale F a , 13 ' 14 from composite models with the composite scale 
F a , 15 from spontaneously broken renormalizable field theories. In the last 
case, the global symmetry must have the gluon anomaly and is called the 
PQ symmetry. 2 If this PQ symmetry is spontaneously broken, there arises 
a pseudo-Goldstone boson 1 coupling to the anomaly with the global sym- 
metry breaking scale F a . 

In QFT, a very light axion is embedded in the phase of a complex SU(2) l 
singlet scalar field s, (it may contain very tiny components (< 10 -7 ) from 
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SU(2) L doublet phases), 4 

s=^J-e ia/F «, a = a + 2%N DW F a , V = N DW F a . (5) 
v2 

So, F a is in general smaller than (s). The potential arising from the anomaly 
term after integrating out the gluon field is the axion potential. Three 
properties of the axion potential are known: 

(i) It is periodic with 2nF a periodicity, 

(ii) The minima are at (a) = 0, 2nF a , 4nF a , ■ ■ ■ , 2 ' 16 

(iii) A set of minima is identical, leaving to a few (Ndw) distinct 
vacua. 17 




The height of the axion potential is the scale A of the nonabelian gauge 
interaction and the boson mixing as shown in Fig. 1. We simply take this 
value as the QCD scale, but in fact it is toAq CD where m is the light quark 
mass. 18 The dominant one Aq CD corresponds to the rf potential. If there 
are quarks, the height is adjusted since as we have seen before a massless 
quark makes it flat. The u and d quark phenomenology gives 

F[a] = Ir ^ F />^(l-cosA). (6 ) 

The essence of the axion solution is that (a) seeks 9 = in the evolving 
universe whatever happened before. It is a cosmological solution 20 as shown 
in Fig. 1. The weak CP violation makes the minimum of the potential 
shifted a little bit at 9 = O(10~ 17 ). The axion mass is given by m a ~ 
(10 7 GeV/ J F Q ) 0.6 eV. 

There are several laboratory experiments, restricting the axion decay 
constants: (i) meson decays, J/^> — * a + 7, T — > aj, K + — > tt + + a, (ii) 
beam dump experiments, p(pr e—)N —> aX,a — > 77 and e + e~, (iii) and 
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nuclear deexcitation, TV* — > Na,a — » 77 and e + e~. Thus, we obtain the 
inner space bound F a > 10 4 GeV from the laboratory experiments. So, 
from the beginning, it was known that the PQWW axion, arising from the 
electroweak scale, is ruled out. 21 Thus, F a has to be very large, having led 
to the so-called invisible axion. But, there is a possibility of detecting it, 22 
and hence it should be called a very light axion or sub-meV axion. 



3. Axion window to outer space 

But the stringent lower bounds on the axion decay constant comes from 
the outer space observations. Firstly, stellar evolutions, if axion existed, 
are affected by axion emissions and the successful standard energy loss 
mechanism due to weak interactions restricts the axion mass toward a 
smaller region, or the axion decay constant to a larger region. The strin- 
gent bound comes from the study of supernova evolution, 23 especially from 
the SN1987A study to give F a > 10 9 GeV. 24 On the other hand, the very 
interesting upper bound on F a is obtained from the axionic contribution to 
dark energy in universe. 20 



3.1. Stars 

The current supernova 24 (globular cluster 25 ) limit on F a is 10 9 GeV(10 10 
GeV) . It uses primarily the Primakoff process with the following coupling, 3 

2 

a 6 ~ 

C ail = - c ajl — -^^F eni F eni => E • B interaction (7) 

C1377 — C-ajj H~ 6 ^ ^ ^iQcm.i — ^0,77 — 1-93, Z = q (8) 
i— light q 

c a77 = determined from high energy physics (9) 

(10) 



l + Z ' ® d — 1 + z 



where the chiral symmetry breaking of u,d quarks are taken into account. 
The number 1.93 corresponds to Z = |. Since the instanton contribution 
to light quark masses is present, 9 ' 10 we may take a band around 1.93. 

In the hot plasma in stars, once produced, they most probably escape 
the core of the star and take out energy. This contributes to the energy 
loss mechanism of star and should not dominate the luminocity: (i) The 
Primakoff process: 7 — > a (present in any model): g ail < 0.6 x 10~ 10 GeV -1 
or F a > 10 7 GeV, and 0.4 eV < m a < 200 kcV ruled out because too heavy 
to produce, 
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io" 5 iff 1 icf 3 icf 2 10 _1 1 10 



Fig. 2. Experimental bound from various experiments, especially including CAST 2003 
data. This figure is taken from the result paper of CAST experiment. 26 Here, we show 
the prediction from Z12 model also. 

(ii) Compton-like scattering: 7e — > ae (DFSZ axion has aee coupling) 
g aee < 2.5 x 1(T 13 , 0.01 eV < m a < 200 keV, and 

(iii) SN1987A, NN -> NNa3 x IO" 10 < g aNN < 3 x IO" 7 =S> F a > 
0.6 x 10 9 GeV. 

Stellar evolution uses the energy loss mechanism, with the aforemen- 
tioned lower bound on F a . But laboratory experiments can offer a more 
effective bound on F a than just the energy loss mechanism, as done by the 
CAST(CERN axion solar telescope) experiment of Fig. 2. 26 

3.2. Universe 

In the standard big bang cosmology (SBB), there is a severe domain wall 
problem. 27 The SBB allows only the domain wall number Now = l- 28 
But the most interesting inflationary cosmology solves this domain wall 
problem at one stroke if the reheating temperature after inflation is lower 
than F a . The inflationary cosmology seems to get support from the COBE 
and WMAP observations of density perturbations in the early universe, 
and NDW problem is not an issue since in SUSY models the reheating 
temperature is required to be smaller than 10 9 GeV. 29 In axion cosmology, 
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the following items are important: 

• The axion decay constant F a , 

• Axion couplings to 7, e,p, n, 

• The domain wall number Npw- 

If a singlet scalar VEV (s) breaks the PQ symmetry, then F a = (s)/Nnw 
defines the domain wall number. Reheating temperature after inflation is 
required to be below F a if Now > 1- 

Axions are created at T ~ F a , but the universe does not change (a) 
until H ~ m a (T = lGeV). Then, the classical field (a) starts to oscillate. 
From the harmonic oscillator type energy density m 2 F 2 , we have m a x 
number density => CDM-like energy: 30 

Pa (T 7 ) = m a (T 7 K(T 7 ) ~ 2.5^^T 3 (^ff (H) 

where the oscillation start-up temperature T\ is the strong interaction scale 
1 GeV and T 7 is the present temperature. If F a is large(> 10 12 GeV), then 
the axion energy density dominates the energy density of the universe. Since 
the energy density is proportional to the number density, it behaves like a 
CDM. 

3.3. Axion window and search for cosmic axions 

The above astro- and cosmological-bounds on F a are summarized as 

10 9 GcV <F a < 10 12 GeV (12) 

If axions are the CDM component of the universe, then they can be 
detected. The feeble coupling can be compensated by a huge number of 
axions. The number density ~ F 2 , and the cross section ~ 1/-F 2 , and there 
is a hope of detecting it. Sikivies cavity detector 22 with dimension of tens 
of cm has been used to give coarse bounds 31 on axion parameters in the 
axion mass of order 10~ 5 eV. 

4. Axions from superstring 

Superstring tells us definite things about global symmetries. If axion is 
present, it is better to be realized in superstring. They are the bosonic de- 
grees in Bmn (MI- axion 13 is B^ v and MD- axion 14 is J3y) and furthermore 
additional massless bosons from compactification are candidates. Super- 
string does not allow global symmetries. But there is an important excep- 
tion to this claim: the shift symmetry of H^ p , which gives the Ml-axion. 
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It is the only allowed global symmetry. Bij are generally heavy; 32 but it is 
a model dependent statement. 

The superstring axion decay constants are expected near the string scale 
which is too large: 33 F a > 10 16 GeV. A key question in superstring models 
is "How can one obtain a low value of F a T' An idea is the following: 

In some compactifications, anomalous U(l) results, 34 where U(l) 
gauge boson eats the Ml-axion to become heavy. Earlier, this di- 
rection, even before discovering anomalous U(l) gauge boson, was 
pointed out by Barr. 35 It became a consistent theory after discov- 
ering the anomalous U(l). Then, a global symmetry survives down 
the string scale. F a may be put in the axion window. It was stressed 
in several references. 36 ' 37 

However, this idea does not work necessarily, as will be commented later. 

Somehow MD-axion(s) may not develop a large superpotential terms. 
But the problem here is the magnitude of the decay constant. MD-axion 
decay constants were tried to be lowered by localizing them at fixed 
points. 19,38 It uses the flux compactification idea and it is possible to have 
a small F a compared to the string scale as in the RS model. One needs the 
so-called throat as schematically shown in Fig. 3. 



Fig. 3. A schematic view of the S2 X S2 X Si throat. At the tip, one elongated S2 can 
shrink to a point shown as a star. The un-shrunk 52 has the cycle C2 where a MD-axion 
resides as a harmonic 2-form. 

Axion mixing Even if we lowered some F a , we must consider hidden sector 
also. In this case, axion mixing must be considered. There is an important 
theorem. 
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Cross theorem on potential heights and decay constants: 39 Suppose 
two axions a\ with F\ and a 2 with F 2 (Fi <C F 2 ) couples to two 
nonabelian groups whose scales having a hierarchy, Ai < A 2 . The 
higher potential A 2 couples to both axions. Then, the diagonaliza- 
tion process chooses the larger potential A 2 corresponds to the 
smaller decay constant and the smaller potential Ai corre- 
sponds to the larger decay constant F 2 . 

So, just obtaining a small decay constant is not enough. Hidden sector may 
steal the smaller decay constant. It is likely that the QCD axion chooses 
the larger decay constant. Recently, the mixing effect has been stressed by 
I.-W. Kim et. al. 19 ' 40 

And most probably, our axion will couple to the \x term: 

H u Hdf(Si, S 2 , ■ ■ ■ )• 

After all, the topologically attractive Bmn may not be the axion we 
want. Let us go back to earlier field theoretic very light axion. In string 
models, its effect toward phenomenology was not calculated before. Now we 
have an explicit model for MSSM, 41 and we can see here whether the idea 
of approximate global symmetry is realized. It must be that at sufficiently 
higher orders the PQ symmetry is broken. In the Z\ 2 _i model, we calculate 
the axion-photon-photon coupling whose result is shown in Fig. 2. But 
the decay constant is at the GUT scale. In this kind of calculation, there 
are so many Yukawa couplings to consider. For example, we encountered 
0(104) terms for d=7 superpotential terms and it is not a trivial task to 
find an approximate PQ symmetry direction. 

In addition, we point out that the MI- axion with anomalous U(l) al- 
ways has a large decay constant since most of the fields are charged under 
this anomalous U(l). Phenomenologically successful axion must need an 
approximate PQ symmetry. 

An approximate PQ global symmetry with discrete symmetry in 
SUGRA was pointed out long time ago: given by Lazarides and Shafi 42 
for a discrete Z 3 x Z 3 . But this field theoretic method does not guarantee 
that string models realize this idea. In this sense, an explicit demonstration 
of an approximate PQ symmetry is vital for a string matter axion. 

5. SUSY extension and axino 

The SUSY extension always introduces gravitino. Gravitinos produced 
thermally after inflation decay very late in cosmic time scale (> 10 3 s) 
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and can dissociate the light nuclei by its decay products. Not to have 
too many gravitinos, the reheating temperature must be bounded, 29,43 
T r < 10 9 GeV (old value), T r < 10 7 GeV (new but model dependent value). 
Therefore, in SUSY theories we must consider the relatively small reheat- 
ing temperature, and the domain wall problem does not matter in axion 
cosmology. 

The SUSY extension with the strong CP solution via axion introduces 
its superpartner axino. Its cosmological significance is in that it can serve 
as keV range warm dark 44 matter or GeV range cold dark matter. 45 Let us 
comment on its CDM possibility. 

For axino to be CDM, it must be stable or practically stable. Without 
the i?-parity conservation, this can not happen. Thus, we require the prac- 
tical i?-parity conservation for the possibility of axino CDM. In addition, 
for axino to be LSP it must be lighter than the lightest neutralino whose 
mass is expected to be around 100 GeV. Thus, the estimation of the axino 
mass is of prime importance. The conclusion is that there is no theoretical 
upper bound on the axino mass and axino mass can be easily in the GeV 
range. Since axion is almost massless, one expects that its superpartners are 
masslcss in the first approximation. Its scalar partner, saxion, obtains the 
mass of order the soft terms after SUSY breaking. Its cosmological effect 
is relatively late decaying nature, adding more photons after its decay. 46 
Regarding mass, saxion is like the other SM SUSY scalars. 

However, the axino mass is intimately related to the SUSY breaking 
scenario and symmetries of the superpotential. The PQ symmetry allows 
the following superpotential 47 

W = fZ(SiS2 - Ft), Z, S U S 2 : singlets 

There also exist SUSY breaking soft terms. Thus, the following potential is 
obtained 

V = |/| 2 (|Si| 2 + |S 2 | 2 )|Z| 2 + (A 1 fS 1 S 2 Z - A 2 fF 2 a Z + h.c.) 

which determines the VEV of Z. Since S\ and 5*2 are of order F a , (Z) is of 
order the A term. Thus, the fermion partners have the mass matrix of the 
form 

m a fF a \ 

m a fF a , m a = f(Z) 

JFafFa / 

The lightest eigenvalue is the axino mass. The others are of order F a . As 
shown above, the axino mass is basically a free parameter, and expected to 
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be smaller than the naive SUSY breaking scale due to the small coupling. 
Most probably, it is lighter than neutralino \. But its mass can be much 
smaller than the SUSY scale as shown from a superpotential of the form, 48 

W = fZ{S!S 2 - X 2 ) + -(X - M f 

o 

where X carries the vanishing PQ charge. This potential is much more 
complicated to analyze. We show ma — 0(A — 2B + C) + 0{m\, 2 / F a ). 
For the standard pattern of soft terms, we have B = A — m 3 / 2 and C — 
A — 2m 3 / 2 - 49 Thus, the axino mass is of order keV. Even the tree level 
axino mass needs the knowledge on the full superpotential, we treat the 
axino as the LSP which is the most probable choice. Its mass is left as a 
free parameter. KeV axinos can be warm dark matter which is thermal relic. 
GeV axinos can be CDM. In this case, the reheating temperature must be 
low. 45 

If gravitino is the next lightest LSP (NLSP), m a < m 3 / 2 < m x , the 
gravitino problem can be resolved, 50 since the thermally produced graviti- 
nos would decay to axino and axion which do not affect BBN produced 
light elements. 

Most probably, \ would be the NLSP, and the thermal production mech- 
anism restricts the reheating temperature after inflation as summarized in 
Fig. 4. At high reheating temperature, thermal production contributes dom- 
inantly in the axino production. Even though the reheating temperature is 
below critical energy density line, there still exists the CDM possibility by 
the non-thermal production (NTP) axinos. Covi et. al. shows 45 

NTP : Flak 2 = ^-Q x h 2 for ma < m x < m 3 / 2 

In Fig. , NTP axinos can be CDM for relatively low reheating temperature 
< 10 TeV, in the region 

10 MeV < m a < m x , NTP axino as CDM possibility. 

The shaded region corresponds to the MSSM models with Q x h 2 < 10 4 , 
but a small axino mass renders the possibility of axino closing the universe 
or just 30% of the energy density. If all SUSY mass parameters are below 
1 TeV, then Q x h 2 < 100 and sufficient axino energy density requires 

{Low reheating is good in view of 
the recent gravitino problem. 
But not good with leptogenesis. 
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Fig. 4. The solid line gives the upper bound from thermal production on the reheating 
temperature as a function of the axino mass. The dark region is the region where non- 
thermal production can give cosmologically interesting results (f2g TP /i 2 ~ 1). 



6. Conclusion 

I reviewed strong CP and axion. In particular, 

• Solutions of the strong CP problem: Nelson-Barr, m u = 0, axion. 
Axion a is the most attractive and plausible solution. 

• Axions can contribute to CDM. Maybe solar axions are easier to 
detect. Most exciting is, it confirms instanton physics by observa- 
tion. 

• Tried to present a superstring matter axion coupling for the first 
time. A QCD axion from superstring may be a window to string. 

• With SUSY extension, O(GeV) axino can be CDM. It is difficult 
to detect this axino from the DM search, but possible to detect at 
LHC as missing energy. 
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